www.gusucode.com > 小波分析理论与Matlab 7实现 > 小波分析理论与Matlab 7实现/code/Ch9/exam9_5.m
% 装载尺度滤波器 load db8; w = db8; subplot(421); stem(w); title('原始尺度滤波器'); % 计算4个滤波器 [Lo_D,Hi_D,Lo_R,Hi_R] = orthfilt(w); subplot(423); stem(Lo_D); title('分解低通滤波器'); subplot(424); stem(Hi_D); title('分解高通滤波器'); subplot(425); stem(Lo_R); title('重构低通滤波器'); subplot(426); stem(Hi_R); title('重构高通滤波器'); % 检验正交性 df = [Lo_D;Hi_D]; rf = [Lo_R;Hi_R]; id1 = df*df' id2 = rf*rf' % 检验二元转换后的正交性 df = [Lo_D 0 0;Hi_D 0 0]; dft = [0 0 Lo_D; 0 0 Hi_D]; zer = df*dft' % 高频和低频解释 fftld = fft(Lo_D); ffthd = fft(Hi_D); freq = [1:length(Lo_D)]/length(Lo_D); subplot(427); plot(freq,abs(fftld)); title('转换模量:低通'); subplot(428); plot(freq,abs(ffthd)); title('转换模量:高通')